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Abstract:  Military  maneuvers  damage  vegetation  and 
compact  and  rut  soils  on  training  lands,  thereby  in¬ 
creasing  the  likelihood  of  hillslope  runoff  and  soil  ero¬ 
sion.  Soil  Freeze-Thaw  (FT)  processes  can  change 
the  hydraulic  geometry  and  roughness  of  vehicular 
ruts  and  reduce  soil  compaction,  which  often  partially 
restores  the  water  infiltration  rate  that  existed  before 
compaction.  The  efficiency  of  these  FT-induced 
'repairs'  depends  on  soil  water  content  and  FT  inten¬ 
sity.  Initial  tests  showed  that  1)  an  experimental  soil 


bin  designed  and  constructed  for  rut  experiments 
allows  acceptable  simulation  of  field  soil  FT,  and  2) 
the  hydraulic  geometry  of  a  rectangular  rill  in  a  fine  silt 
soil  with  an  initial  volumetric  water  content  of  36% 
changes  dramatically  due  to  rill  sideslope  slumping 
during  thaw.  Future  experiments  will  compare  differ¬ 
ences  in  the  response  of  natural  rills  and  vehicular  ruts 
to  FT-induced  soil  failure,  and  investigate  the  effects  of 
FT  on  soil  erodibility  and  the  influences  of  snow  cover 
on  soil  erosion  processes  in  the  spring. 
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NOMENCLATURE 

A  =  cross-sectional  flow  area 
C,n-  roughness  coefficients 
d  =  rill  cross-sectional  depth 
F  =  freeze 
FI  =  first  freeze 
F2  =  second  freeze 

FERF  =  Frost  Effects  Research  Facility 
FT  =  freeze-thaw 
FTC  =  freeze-thaw  cycle 
Q  =  runoff  volume 
R  =  hydraulic  radius 
RG  =  resistivity  gage 


S  =  channel  bottom  slope 
T  =  thaw 
T1  =  first  thaw 
T2  =  second  thaw 

Ta  =  air  temperature  between  freezing  panels 
and  soil  surface  in  bin 

Tp  =  panel  temperature 

Tssc  =  temperature  of  the  soil  next  to  RG  C 

Tssr  =  temperature  of  the  soil  in  mid-rill 

V  =  flow  velocity 

xv  =  rill  width 
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INTRODUCTION 
Soil  erosion 

Soil  is  naturally  eroded  by  water  flowing 
down  bare  or  partially  vegetated  hillslopes;  this 
erosion  is  a  function  of  the  erodibility  (detachabil- 
ity)  of  soil  particles  and  the  transport  power 
of  flowing  water  (Table  1).  Soil  erodibility  is  a 
function  of  interparticle  friction,  bonding,  and 
interlocking-  The  strength  of  these  particle  inter¬ 
actions  depends  on  soil  particle  size  and  distribu¬ 
tion,  soil  structure  and  structural  stability,  soil 
permeability,  water  content,  organic  matter  con¬ 
tent,  and  clay,  mineral,  and  chemical  constituents 
(Lai  and  Elliot  1994).  Crusts  that  often  form  on 
soil  surfaces  also  increase  the  resistance  of  soil 
particles  and  aggregates  to  erosion. 

In  addition,  vegetation  cover  significantly  in¬ 
fluences  the  amount  of  soil  erosion  at  a  location 
by  1)  binding  surface  and  near-surface  soil  parti¬ 


cles  through  root-soil  bonds,  2)  reducing  rainfall 
erosivity  by  intercepting  free-falling  raindrops 
and  reducing  the  number  and  kinetic  energy  of 
raindrop  impacts  onto  the  soil  surface  (Evans 
1980),  and  3)  reducing  overland  flow  velocities  by 
increasing  friction  where  plant  stems  protrude 
into  the  flow  (Prosser  et  al.  1995).  Thus,  any  pro¬ 
cess  that  alters  either  the  physical  strength  of  a 
soil  or  reduces  its  vegetative  cover  increases  the 
likelihood  of  water  erosion  of  that  soil. 

Effects  of  Army  maneuvers  on 
training  land  erosion 

The  Army  is  responsible  for  over  18,500  square 
miles  (4.8  x  104  km2)  of  training  lands  that  are 
principally  used  to  ensure  the  military  readiness 
of  its  units  (Doe  1992).  Trainers  must  maintain 
realism  during  maneuvers,  but  the  goals  of  readi¬ 
ness  and  realism  often  conflict  with  environmen- 


Table  1.  Factors  that  determine  the  severity  of  water  erosion.  (After  Lai  1994.) 

I.  Climatic  erosivity 

1.  Rainfall  erosivity 

•  A  measure  of  the  ability  of  rain  to  detach  sediment  particles  and  surpass  the  infiltration  capacity  of  soil  so  that 
overland  flow  begins. 

•  Also  called  the  Energy-Intensity  (El)  parameter,  a  function  of  rainfall  volume,  raindrop  impact,  and  peak  intensity. 
2.  Runoff  erosivity 

•  A  measure  of  the  ability  of  flowing  water  to  detach  and  transport  sediment  particles. 

•  A  function  of  runoff  volume  and  peak  flow. 

II.  Soil  erodibility 

•  A  measure  of  the  susceptibility  of  sediment  particles  to  being  detached  and  transported  by  rain  and  flowing  water. 

•  A  function  of  soil  texture,  structure,  permeability,  organic  matter  content,  chemical  constituents,  and  clay 
mineralogy. 

III.  Topography 

•  Hillslope  length,  steepness,  and  shape  influence  overland  flow  velocities  and  turbulence,  which  partially  deter¬ 
mines  the  likelihood  of  rill  formation. 

IV.  Land  use 

•  Disturbance  to  a  soil  surface  influences  the  effectiveness  of  raindrop  impacts  in  moving  soil  particles,  soil  infiltration 
rates,  and  overland  flow  velocities  and  turbulence. 
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Figure  1.  Runoff  power  vs.  sediment  yield  from 
1-m2  erosion  plots  used  and  not  used  by  off-road 
vehicles .  (After  Iverson  1980.) 

tal  requirements.  Army  land  managers,  who  need 
to  preserve  natural  resources  on  training  lands,  are 
required  to  minimize  vegetation  damage  and  soil 
disruption,  which  inevitably  occur  during  train¬ 
ing.  Vegetation  and  soils  are  important  resources 
in  themselves,  but  damage  to  them  often  leads  to 
accelerated  soil  erosion  (Dregne  1983). 

An  armored  unit  on  the  move  or  an  infantry 
march  damages  vegetation,  breaks  up  soil  crusts, 
loosens  surface  soil,  alters  soil  structure,  weakens 
soil  aggregates,  changes  soil  surface  roughness, 
changes  the  shape  and  number  of  surface  depres¬ 
sions,  and  often  compacts  soils.  Compaction  in¬ 
creases  soil  bulk  density,  reduces  infiltration  and 
hydraulic  conductivity  (Voorhees  et  al.  1979;  Webb 
1983;  Thornes  1980;  Braunack  1986a, b;  Ayers  1994; 
Campbell  1994;  Horton  et  al.  1994),  and  restricts 
soil  aeration,  which  impairs  root  growth,  plant 
nutrient  uptake,  and  seedling  emergence  (Chan¬ 
cellor  1977,  van  Ouwerkerk  1991,  Stepniewski  et 
al.  1994).  Thurow  et  al.  (1993)  report  that  recovery 
of  damaged  plants  and  new  growth  is  often  limit¬ 
ed  on  compacted  soils.  This  results  in  a  vegetative 
cover  that  is  too  sparse  or  composed  of  species  less 
effective  in  protecting  and  binding  soil  particles 
sufficiently  to  contribute  to  their  stability  and  flow 
resistance. 


The  impact  of  soil  compaction  on  hillslope  hy¬ 
drology  and  erosion  processes  is  substantial. 
Rainfall  of  a  given  intensity,  which  usually  infil¬ 
trates  into  undisturbed  soils,  often  does  not  infil¬ 
trate  into  the  same  soils  if  compacted.  Vehicular 
ruts  formed  during  maneuvers  often  have  com¬ 
pacted  soils  beneath  and  adjacent  to  them.  This 
additional  surface  water  adds  to  runoff  volumes 
(Eckert  et  al.  1979,  Mathier  and  Roy  1993)  and 
makes  runoff  periods  longer  (Hinckley  et  al. 
1983).  I  believe  that  the  ruts  tend  to  channel  the 
additional  runoff  into  rill-type  flows,  with  veloc¬ 
ities  that  may  be  measurably  higher  than  veloc¬ 
ities  in  natural  rills  on  the  same  hillslope.  Natural 
rills  don't  carry  any  additional  runoff  coming 
from  compacted  soils.  I  suggest  that  such  higher 
runoff  erosivity  in  ruts  may  explain  why  gullies 
on  training  lands  can  form  and  enlarge  faster 
than  they  do  on  adjacent  undisturbed  soils. 

Iverson  (1980)  showed  that,  for  a  given  runoff 
power,  more  sediment  was  eroded  from  hillslope 
plots  that  were  used  by  off-road  vehicles  than 
from  those  that  were  not  used  (Fig.  1).  He  report¬ 
ed  that  the  flow  capacity  on  used  plots  increased 
more  than  linearly  with  runoff  power  because  of 
increased  runoff  volume  and  flow  channeliza¬ 
tion.  Reduced  infiltration  and  frictional  resistance 
to  flow  on  used  plots  cause  overland  runoff  to 
happen  more  rapidly  and  attain  an  eroding  dis¬ 
charge  over  a  larger  portion  of  a  used  hillslope 
than  on  an  unused  slope. 

EFFECTS  OF  GROUND  FREEZING 
ON  MANEUVER  IMPACTS 

Soil  compaction 

Soil  compaction  is  the  compression  of  unsatur¬ 
ated  soil  because  of  reduction  of  its  air-filled  pore 
space  without  a  change  in  mass  wetness.  Com¬ 
paction  results  from  simultaneous  application  of 
vertical  pressures  and  shearing  stresses  from  traf¬ 
ficking  on  soils  (Hillel  1980).  The  amount  of 
vehicular  soil  compaction  is  determined  by  vehi¬ 
cle  type,  pattern  of  loading  (static,  dynamic,  sta¬ 
ble,  vibratory),  vehicle-traffic  motion  (straight  or 
turning),  number  of  vehicle  passes,  vehicle  pres¬ 
sures  applied,  soil  texture,  density  and  moisture, 
and  state  and  stability  of  soil  structure  (Voorhees 
et  al.  1978,  1986;  Akram  and  Kemper  1979;  Hillel 
1980;  Webb  1983;  Gupta  et  al.  1989;  Braunack 
1986a,b;  Foltz  1992;  Thurow  et  al.  1993). 

Depending  on  the  interplay  of  these  factors,  a 
compacting  vehicular  force  (applied  load) 
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Figure  2.  Water  retention  characteristics  of  Nicollet  clay  loam  at  vari¬ 
ous  bulk  densities;  dense  soils  hold  less  xuater  at  high  matric  potential 
and  more  ivater  at  low  potential.  (After  Gupta  et  al.  1989.) 


•  Establishes  or  breaks  bonds  in  aggregates 
and  among  particles  (Koolen  and  Kuipers 
1983). 

•  Crushes  soil  aggregates,  thus  making  the  soil 
particle  distribution  more  uniform,  which 
tends  to  reduce  soil  pore  size  distribution, 
which  changes  the  soil  water  suction  and 
water  retention  characteristics  (Taylor  and 
Box  1961,  Koolen  and  Kuipers  1983,  Gupta  et 
al.  1989,  Cruse  and  Gupta  1991)  (Fig.  2). 

•  Increases  the  volumetric  water  content  be¬ 
cause  the  applied  stresses  act  for  a  very  short 
time  during  compaction,  which  strongly 
restricts  the  amount  of  water  that  moves  out 
of  the  soil,  i.e.,  traps  the  water  (Hillel  1980, 
Koolen  and  Kuipers  1983). 

•  Reduces  the  volume  of  soil  voids,  which 
reduces  total  porosity  and  increases  dry  bulk 
density  (Hillel  1980,  Kooistra  and  Tovey 
1994). 

•  Changes  soil  pore  geometry  (Gupta  et  al. 
1989)  and  reduces  the  interconnectedness  of 
larger  pores  (Hillel  1980,  Iverson  1980). 

•  Increases  soil  penetration  resistance  (Voor- 
hees  et  al.  1986),  especially  with  rubber  tires, 
not  tracks  (Braunack  1986b). 

•  Decreases  soil  infiltration  (Akram  and  Kem¬ 
per  1979)  (Fig.  3)  and  permeability  (Braunack 
1986a,b);  and  if,  during  trafficking,  the  vehi¬ 
cles'  wheels  slip,  realigns  soil  particles  paral¬ 
lel  to  the  direction  of  the  shear  forces,  which 
causes  additional  compaction  and  further  re¬ 
duces  infiltration  (Gupta  et  al.  1989). 


•  Increases  surface  runoff  and  accelerates  soil 
erosion  (Foltz  1992,  van  Ouwerkerk  and 
Soane  1994). 

Compacted  soil  particles  can  be  loosened  by  the 
shrinking  and  swelling  associated  with  wetting 
and  drying  of  clays  (Larson  and  Allmaras  1971), 
by  root  growth  and  by  soil  freeze-thaw  (FT),  so 
that  the  infiltration  rate  of  a  previously  compacted 
soil  often  returns  to  or  near  to  its  previous  values 
(Canarache  1991,  Thurow  et  al.  1993).  The  ice  that 
forms  in  soil  pores  during  freezing  can  reduce  the 
density  of  compacted  soil  by  pushing  soil  grains 
apart  and  reducing  their  degree  of  interlocking. 
This  is  most  effective  when  volumetric  soil  water 
content  approaches  soil  porosity.  The  amount  of 
soil  ice  formed  in  pores  or  as  lenses  is  related  to  the 
volume  of  soil  water  present  when  freezing  starts 
and  the  volume  drawn  to  the  freezing  zone  from 
the  subsoil  (Miller  1980). 

Upon  thawing,  a  soil  is  usually  less  dense, 
although  the  degree  of  FT-induced  soil  expansion 
depends  on  soil  water  content,  soil  texture  and 
depth,  the  rate  of  frost  penetration,  the  number  of 
FT  cycles,  and  the  depth  of  compaction  (Webb  et 
al.  1983).  This  soil  loosening  can  be  sufficient  to 
increase  infiltration  and  reduce  runoff,  especially 
in  arid  and  semi-arid  environments  (Schumm  and 
Lusby  1963).  However,  research  results  to  date 
conflict  regarding  the  degree  of  FT-induced  loos¬ 
ening,  as  discussed  below. 

Chamberlain  and  Gow  (1979)  found  that  FT 
reduced  void  ratio  and  increased  vertical  perme¬ 
ability,  while  Blake  et  al.  (1976)  found  that  FT  did 
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Figure  3.  Effects  of  compacting  loads  and  water  content  on  infiltration  rate,  compaction,  and  bulk  density.  ( After 
Akram  and  Kemper  1979.) 


not  reduce  bulk  density  of  compacted  Nicollet 
clay  loam  in  Minnesota  at  20-30  cm  depth,  even 
after  9  years.  Akram  and  Kemper  (1979)  reported 
that  FT  reduced  soil  compaction  in  different  soils 
and,  as  a  result,  infiltration  rates  increased  after 
FT.  Infiltration  increased  most  after  one  FT  cycle  in 
a  loamy  sand,  after  two  cycles  in  one  sandy  loam 
and  three  cycles  in  another  sandy  loam,  and  infil¬ 
tration  was  still  increasing  in  a  clay  loam  after  four 
cycles  (Fig.  4). 

Voorhees  et  al.  (1978,  1986)  suggested  that  the 
persistence  of  soil  compaction  in  agricultural 
fields,  despite  FT  cycling,  may  be  partly  attribut¬ 
able  to  increasing  tractor  weight.  FT  loosened  the 
soil  to  about  20  cm  but  compaction  persisted  be¬ 
low  20  to  90  cm  depth  where  it  was  unaffected  by 
FT  (Fig.  5);  however,  this  persistence  was  in  part 
soil-dependent.  Voorhees  (1983)  reported  that  soil 
FT  and  wetting  and  drying  reduced  soil  pene¬ 
tration  resistance  by  20-50%  and  that  the  FT  was 
more  effective  in  reducing  resistance  when  the  soil 
was  wetter  at  freezeup. 


The  above  results  may  differ  because  of  the  real 
variability  in  nature  and  the  absence  of  standard¬ 
ized  methods  for  characterizing  soil  compactness 
(Soane  and  van  Ouwerkerk  1994),  although  soil 
bulk  density,  total  porosity,  void  ratio,  specific  vol¬ 
ume,  and  unit  weight  are  generally  considered  to 
be  fundamental  criteria  that  define  the  degree  of 
compaction. 

Vehicular  ruts 

Vehicles  can  rut  a  soil  surface,  depending  on 
vehicle  load  in  relation  to  soil  conditions  at  the  time 
of  trafficking  (Richmond  et  al.  1995).  Ruts  are 
potential  sites  of  high  soil  erosion  when  aligned 
directly  up  and  down  slope  and  are  hydraulically 
similar  to  natural  rills  in  the  erosivity  of  the  flows 
in  them.  Water  flow  is  faster  and  more  turbulent  in 
natural  rills  than  it  is  in  overland  sheet  wash 
(Evans  1980)  and  it  has  more  energy  to  detach  and 
transport  sediment. 

Voorhees  et  al.  (1979)  reported  that  the  ruts  left 
by  wheeled  vehicles  can  act  as  channels  to  concen- 
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Following  Compaction 

Figure  4.  Effect  of  FT  cycles  on  infiltration  rate  of 
compacted  soils.  (After  Akram  and  Kemper  1979.) 

trate  surface  runoff,  which  increases  its  sediment 
transport  capacity  and  subsequent  soil  erosion. 
Foltz  (1993)  determined  that  there  is  200-400% 
more  erosion  on  rutted  roads  than  on  unrutted 
roads.  Morgan  (1977)  determined  that  rill  sediment 
transport  exceeded  that  on  inter-rill  surfaces  by  a 
factor  of  40  on  an  11°  slope.  Mutchler  and  Young 
(1975)  determined  that  more  than  80%  of  eroded 


hillslope  sediment  is  transported  in  rills,  and 
Meyer  et  al.  (1975)  reported  a  threefold  increase 
in  soil  loss  following  rill  development  on  a  hill- 
slope.  Thus,  rill  (and  rut)  erosion  is  a  geomorphi- 
cally  significant  process  (Slattery  and  Bryan 
1992). 

The  intermittent  flows  in  rills  and  ruts  are  usu¬ 
ally  appreciably  deeper  than  the  height  of  the 
coarsest  roughness  elements  within  them  and,  in 
this  regard,  are  hydraulically  similar  to  conven¬ 
tional  open-channel  flows  in  rivers  (Thornes 
1980,  Hairsine  and  Rose  1992).  The  erosivity  of 
such  flows  is  directly  related  to  runoff  volume 
and  velocity 

where  Q  =  runoff  volume 

A  =  cross-sectional  flow  area 
V  =  flow  velocity. 

The  velocity  is  determined  by  channel  roughness, 
cross-sectional  shape,  and  slope  and  can  be  esti¬ 
mated  by  the  Chezy  (eq  1)  or  Manning's  (eq  2) 
equations 

V  =  CjRS  (1) 

V=  1.49/tt  R2/3  S1'2  (2) 

where  V  =  flow  velocity 

C,  n  =  roughness  coefficients 

R  =  hydraulic  radius,  A/P  (P,  wetted  peri¬ 
meter) 

S  =  channel  bottom  slope  (approximately). 

It  is  apparent  that  processes  that  change  rill  or 
rut  cross-sectional  shape  can  have  a  major  influ¬ 
ence  on  flow  erosivity  and  thus  on  sediment 


Figure  5.  Penetrometer  resistance  4  years  after  initial  loading  of  Webster  clay 
loam  in  Waseca,  Minnesota ,  with  a  wheeled  vehicle.  (After  Voorhees  et  al  1986.) 
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Figure  6.  Hypothetical  changes  in  rill  hydraulic  geometry 
because  ofFT-induced  soil  failures;  as  thawed  soil  along  rill 
side  slopes  collapses ,  the  rill  cross-sectional  depth  (d) 
decreases  and  width  (w)  increases ,  and  the  character  of  the 
rill  flow  becomes  more  like  that  of  less-erosive  overland  flow. 
This  evolution  assumes  no  intermittent  flow  to  clean  out  the  rill 
channel 


Figure  7.  Sediment  concentrations  in  runoff  over  variably  compacted 
soil  was  higher  in  denser  soil  because  the  surface  roughness  of  and  the 
infiltration  into  the  more  compacted  soil  was  less ,  resulting  in  high , 
near-bed  water  velocities  and  more  bed  soil-particle  detachment  and 
transport.  (After  Parker  et  al.  1995.) 


transport  and  hillslope  erosion.  In  some  loca¬ 
tions,  frost-induced  soil  creep  is  more  effective  at 
moving  surface  soils  downslope  than  is  summer 
creep  and  can  obliterate  natural  rills  over  one 
winter  (Schumm  and  Lusby  1963;  Schumm  1956, 
1964,  1967;  Carson  and  Kirkby  1972).  Such  rills 
become  clogged  with  sediment  and  the  initial 
runoff  in  the  spring  often  becomes  sediment 
laden  as  it  clears  out  infilled  rills  (Schumm  and 
Lusby  1963).  Thus,  rill  enlargement  often  results 


from  sediment  infilling  processes  and  sediment 
clearing  and  channel  cutting  by  rill  flows  (Kirkby 
1980,  Piest  et  al.  1975).  Owoputi  and  Stolte  (1995) 
state  that  further  study  of  the  characteristics  of 
rills  is  necessary  for  an  improved  prediction  of 
the  soil  erosion  process. 

I  hypothesize  that  the  cross-sectional  shape  of 
a  rill  would  become  broader  and  shallower  after 
FT-induced  failures  along  the  rill  side  slopes  and 
no  intermittent  flows  (Fig.  6).  As  a  rill's  depth 


decreases,  bankfull  flow  within  it  tends  toward 
sheet  flow,  where  width  is  much  greater  than 
depth.  Thus,  the  amount  of  hillslope  erosion  tak¬ 
ing  place  in  such  a  rill  could  be  less  than  that  in  a 
rill  unaltered  by  FT.  In  nature,  sediment  infilling 
and  channel  alterations  by  FT-induced  soil  creep 
occur  in  conjunction  with  sediment  transport  by 
intermittent  flows,  and  rill  cross-sectional  shape  is 
determined  by  hydraulic  and  soil  processes  in 
northern  climes. 

I  am  unaware  of  studies  that  have  compared 
erosion  on  a  compacted  and  rutted  hill  to  that  on 
an  undisturbed,  rilled  slope,  but  I  hypothesize  that 
the  transport  capacity  of  vehicular  ruts  may  be 
changed  by  FT  processes  as  are  rills.  However,  the 
ruts  may  not  infill  as  rapidly  as  natural  rills  be¬ 
cause  the  compacted  soil  along  a  rut  may  have 
higher  shear  strength,  which  may  temporarily  re¬ 
tard  FT-induced  creep  down  the  rut  side  slope. 
Thus,  I  also  hypothesize  that  erosion  on  a  rutted 
hill  could  be  more  severe.  First,  flows  in  ruts  may 
have  higher  velocity  and,  thus,  higher  stream 
power  than  rill  flows.  Parker  et  al.  (1995)  found 
that  flow  velocity  near  the  soil  bed  increased  with 
bulk  density  (compaction),  because  surface  rough¬ 
ness  was  less  in  the  more  compacted  soil,  and  that 
sediment  concentration  in  runoff  was  higher  on 
more  highly  compacted  soil  (Fig.  7).  They  found 
that  the  increased  velocity  from  compaction  had 
more  effect  on  erosion  than  the  increased  soil 
shear  strength.  Second,  the  volume  of  surface  wa¬ 
ter  would  be  greater  on  a  compacted  slope  than  on 
the  undisturbed  slope,  owing  to  reduced  infiltra¬ 
tion.  Third,  the  compacted  soil  around  a  rut  would 
allow  less  infiltration  along  its  length  than  may 
occur  along  a  rill  on  the  same  slope;  thus,  the  ero- 
sivity  of  rut  flow  could  be  maintained  for  greater 
distances  downslope.  Research  comparing  rill  and 
rut  processes  is  planned  as  part  of  this  project. 

Seasonal  soil  erodibility 

While  FT  can  loosen  compacted  soil  and 
smooth  rills,  and  possibly  ruts,  over  time,  it  can 
also  make  undisturbed  soils  more  erodible  in  the 
spring  than  they  are  at  other  times  of  the  year 
(Gatto  1995).  When  soil  water  freezes  in  the  win¬ 
ter,  the  ice  crystals  can  disrupt  soil  grain  interlock¬ 
ing,  which  results  in  a  less  dense,  weak  soil  upon 
thaw.  In  addition,  the  soil  water  content  in  a  freez¬ 
ing  soil  usually  increases  as  moisture  is  drawn 
from  the  unfrozen  soil  below  to  the  freezing  front. 
Thus,  newly  thawed  soil  usually  has  more  water 
than  before  it  froze  (it  may  be  temporarily  satur¬ 
ated),  which  contributes  to  its  low  strength  and 


makes  it  highly  susceptible  to  downslope  move¬ 
ment  by  gravity,  detachment  by  raindrops,  and 
detachment  and  transport  by  overland  flow. 

The  magnitude  of  these  FT-induced  effects  is 
variable,  however;  Benoit  and  Voorhees  (1990), 
and  Kok  and  McCool  (1990)  report  that  FT  effects 
are  some  of  the  least  understood  aspects  of  the  soil 
erosion  process,  even  though  soil  FT  processes 
have  been  investigated  for  years.  In  addition,  cur¬ 
rent  erosion  models  still  cannot  predict  infiltrabil- 
ity  or  erodibility  in  frozen  and  partially  frozen 
soils  at  their  surface  (Seyfried  and  Flerchinger 
1994),  even  though  Young  et  al.  (1993)  have  devel¬ 
oped  a  method  to  predict  soil  frost  depth  as  part  of 
the  Water  Erosion  Prediction  Project  (WEPP)  soil 
erosion  model. 

Knowledge  of  the  role  of  snow  cover  on  soil 
erosion  mechanics  remains  rudimentary  as  well. 
Haupt's  (1967)  research  showed  that  a  snow  cover 
tends  to  insulate  and  thus  preserve  soil  frost,  even 
during  spring  rains,  by  preventing  raindrops  from 
contacting  the  frozen  soil.  His  study  shows  that 
because  the  soil  remains  frozen  under  snow- 
covered  plots  and  is  very  resistant  to  flow,  virtu¬ 
ally  no  soil  is  eroded,  even  though  the  vegetation 
below  the  snow  is  sparse. 

A  significant  portion  of  previous  soil  erosion  re¬ 
search  has  been  directed  at  changes  in  soil  bulk 
properties,  not  changes  in  surface  soil  structure 
and  strength,  which  are  most  important  to  erodi¬ 
bility.  Shainberg  et  al.  (1994)  point  out  that  soil 
detachment  by  rill  or  overland  flows  depends  on 
soil-particle-binding  forces  at  the  soil/water  inter¬ 
face,  not  bulk-tensile-strength  properties  of  the 
soil  at  depth.  Yet,  bulk  soil  strength  properties  are 
often  measured  and  used  to  predict  soil  erodibil- 
ity. 

•  Misra  and  Rose  (1995)  summarized  research  on 
the  relationships  between  soil  strength,  as  mea¬ 
sured  in  the  field,  and  soil  erodibility,  as  defined 
by  the  amount  of  soil  particles  moved  by  rain 
splash  and  runoff.  However,  that  relationship 
remains  unclear,  and  the  need  for  a  method  to  pre¬ 
dict  seasonal  soil  erodibility  over  time  persists, 
which  impedes  improvements  in  our  ability  to 
predict  soil  erosion  (Nearing  et  al.  1994). 

RESEARCH  NEEDED  AND  PROJECT  GOALS 

Clearly,  soil  freeze-thaw  cycling  is  a  dynamic 
process,  substantially  affecting  runoff  and  soil 
erodibility  during  the  year.  And  the  most  impor¬ 
tant  experimental  topics  in  the  soil  erosion  arena 
are  the  dynamic  processes  that  determine  soil 
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resistance  to  hydraulic  forces  and  determine  run¬ 
off  formation  rates  and  volumes  (Kirkby  1980, 
Gerits  et  al.  1990).  Papendick  and  Saxton  (1990)  re¬ 
ported  that  research  on  frozen  soil  effects  remains 
a  high  priority.  Cooley  (1990)  described  the  impor¬ 
tance  of  incorporating  the  effects  of  FT  on  soil 
compaction  and  strength,  runoff,  and  erosion  into 
soil  erosion  models.  In  spite  of  the  accelerated  ero¬ 
sion  that  often  results  from  Army  maneuvers,  past 
research  has  not  determined  the  significance  of  FT 
in  alleviating  vehicular  compaction,  in  reestab¬ 
lishing  soil  infiltration,  in  changing  rut  geometry, 
and  in  determining  sediment-transport  capacity 
on  training  lands. 

Available  research  results  conflict  regarding  the 
efficacy  of  FT  in  reducing  vehicular  compaction, 
do  not  define  the  basic  soil-FT  processes  involved, 
and  are  absent  regarding  FT  effects  on  vehicular 
ruts.  The  effects  of  FT  on  natural  rills  has  been  in¬ 
vestigated  at  a  preliminary  level.  And  we  know 
that  snow  cover  can  retard  soil  thaw  and  absorb 
raindrop  impact  energy,  but  we  lack  details  on  how 
snow  affects  soil  particle  detachment  on  the  soil  sur¬ 
face  and  affects  soil  creep  along  rills  and  ruts. 

Laboratory  and  field  experiments  will  be  con¬ 
ducted  to 

•  Measure  FT-induced  changes  in  shear 
strength,  penetration  resistance,  infiltration 
rates,  and  surface  geometry  and  roughness  of 
compacted  and  rutted  soils  for  different  FT 
regimes,  compaction  loads,  and  soil  type  and 
water  content. 

•  Determine  if  compaction  and  ruts  from 
tracked  and  wheeled  vehicles  are  affected 
differently  by  soil  FT. 

•  Determine  differences  in  the  cross-sectional 
shapes  of  rills  and  ruts  caused  by  FT-induced 
soil  creep. 

•  Evaluate  how  and  to  what  degree  soil  FT 
rates  and  soil  water  contents  affect  soil  shear 
strength  and  penetration  resistance  (soil 
erodibility)  of  different  soil  types. 

•  Use  rainfall  simulators  to  determine  soil 
erodibility  with  various  antecedent  soil 
water  conditions  and  soil  types,  and  rainfall 
erosivity  before  frost  and  after  thaw  on  dif¬ 
ferent  slopes. 

•  Determine  if  snow  cover  alters  the  tempera¬ 
ture  gradient  in  the  soil  sufficiently  to  reduce 
the  amount  of  water  drawn  to  a  soil  freezing 
zone. 

•  Evaluate  the  effects  of  snow  cover  on  the 
amount  of  thaw  creep  along  rills  and  ruts. 

•  Investigate  soil  particle  detachment  under 
snow. 


Laboratory  studies  will  include  large-  and 
small-scale  experiments  carried  out  in  the  Frost 
Effects  Research  Facility  (FERF)  at  CRREL  and 
small-scale  experiments  in  coldrooms.  Field 
research  will  be  done  in  three  hydro-climatic 
regions:  the  cool,  semi-arid  zone  of  south-central 
Washington  at  the  Army's  Yakima  Training  Center 
(YTC),  the  cold,  humid  area  of  upper  New  Eng¬ 
land  at  CRREL  or  Ethan  Allen  Firing  Range 
(EAFR),  Vermont,  and  the  cold,  dry  climate  of  the 
upper  midwest,  possibly  at  Camp  Ripley,  Minne¬ 
sota,  or  Fort  McCoy,  Wisconsin.  The  studies  at  YTC 
and  EAFR  are  currently  underway. 

This  research  will  extend  existing  knowledge  of 
soil  erosion  mechanics,  developed  for  the  tradi¬ 
tionally  agricultural  and  rangeland  settings,  to 
military  training  lands  where  soil  frost  forms.  It 
will  provide  information  to  terrain  modelers  to 
help  improve  their  simulations  of  seasonally 
dynamic  soil  processes  that  significantly  influence 
terrain  evolution  and  hillslope  soil  processes.  The 
research  results  may  eventually  be  used  to  modify 
the  Revised  Universal  Soil  Loss  Equation  (RUSLE) 
and  WEPP  soil  erosion  models  to  better  simulate 
the  effects  of  winter  processes  and  conditions  on 
runoff  erosivity  and  soil  erodibility,  information 
that  is  now  unavailable  to  Army  land  managers. 

OBJECTIVES  OF  THE  INITIAL  TESTS 

The  objectives  in  these  tests  were  1)  to  see  if  I 
could  achieve  the  freeze  rates  in  the  soil  bin  locat¬ 
ed  in  the  FERF  that  would  allow  me  to  do  my 
planned  experiments,  2)  to  determine  if  the  soil  in 
the  bin  would  freeze  and  thaw  in  a  natural  fashion 
and,  thus,  reasonably  simulate  field  conditions,  3) 
to  evaluate  the  operation  and  adequacy  of  the  in¬ 
struments  to  be  used  in  future  experiments  with 
multiple  bins,  and  4)  to  measure  cross-sectional 
changes  in  a  simulated,  rectangular  rill. 

APPROACH 
Bin  construction 

The  first  bin  was  constructed  with  two  layers  of 
pressure-treated,  3/4-in.  (1.9-cm)  plywood,  rein¬ 
forced  with  pressure-treated,  2-  x  6-in.  (5-  x  15-cm) 
boards  and  insulated  with  one  layer  of  2-in.  (5-cm) 
polystyrene  foam  board  on  the  sides  and  bottom 
(Fig.  8a).  The  inside  of  the  bin  was  lined  with  an  im¬ 
permeable  membrane  (Fig.  8b)  that  was  trimmed 
and  tacked  to  the  interior  walls.  The  bin  in  final  posi¬ 
tion  for  these  first  tests  was  horizontal  along  its  long 
axis  but  sloped  from  its  east  side  to  its  west  approx¬ 
imately  0.5  in.  in  4  ft  (1.3  cm  in  1.2  m).  The  bin  is  a 
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a.  General  construction. 


b.  Bin  liner  (south  side  in  foreground). 
Figure  8.  Experimental  soil  bin. 


closed  system,  so  the  only  water  in 
the  soil  during  the  tests  is  that 
present  at  the  start. 

Soil 

The  soil  is  a  low-plasticity,  inor¬ 
ganic  clayey-silt,  with  81.6%  silt- 
and  clay-sized  particles  and  17.7% 
fine  sands  (Fig.  9a)  and  is  classified 
as  ML  in  the  Unified  Soil  Classifica¬ 
tion  System.  It  has  a  specific  gravity 
of  2.72,  a  liquid  limit  of  28%,  and  a 
plastic  index  of  1  (Shoop  and  Gatto 
1992).  The  Proctor  test  shows  that  its 
maximum  dry  density  is  achieved  at 
17%  gravimetric  soil  water  content 
(Fig.  9b).  The  silt  was  wetted  to  a 
gravimetric  water  content  of  21  to 
23%  (Table  2),  which  is  92  to  97%  sat¬ 
urated  to  simulate  a  thawed,  satu¬ 
rated  silt  at  a  cold,  humid  site  in  ear¬ 
ly  spring,  and  placed  into  the  bin. 

We  placed  two  to  three  backhoe 
buckets  of  soil  in  the  bin,  spreading 
each  bucket  load  with  a  shovel  and 
raking  it  level.  When  this  loose  soil 
was  about  8  in.  (20  cm)  thick,  we 
tamped  it  to  a  6-in.  (15-cm)  soil  layer 
and  roughened  the  surface  with  a 
rake  before  the  next  backhoe  load 
was  placed  to  reduce  the  develop¬ 
ment  of  boundaries  between  the  lay¬ 
ers.  Three  soil  layers  were  prepared 
in  this  fashion.  Dry  densities  and 
volumetric  water  contents  of  the  soil 
at  the  surface  of  each  layer  were  103- 
104  lb/ft3  (1.64-1.67  Mg/m3)  and 
35-38%,  respectively  (Table  2). 

On  the  surface  of  the  final  soil 
layer,  we  formed  a  rectangular  rill 


COBBLES 

GRAVEL 

SAND 

SILT  or  CLAY 

Crs  1  Fine 

Crs  1  Medium!  Fine 

U.S.  Sieve  Size  (in) 

U.S.  Standard  Sieve  # 

Hydrometer 

a.  Grain  size  distribution  (ASTM  1996a). 


b.  Proctor  test  results  with  5.5-lb  (2.5-kg)  ham¬ 
mer  (ASTM  1996b). 

Figure  9.  Silt  used  in  tests. 
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Percent  Retained  by  Weight 


Table  2.  Initial  conditions  on  the  surface  of  each  soil  layer.  Cores  taken  and 
calculations  made  by  Richard  Roberts,  CRREL,  25  April  1995.  Sample  locations 
shown  in  Figure  13.  Core  size  =  0.01  ft3  =  283.2  cm3. 


Layer/core  no. 

_ Z _ 

(lb/ft3)  (Mx/m3) 

(lb/ft3) 

Td _ 

(Mg/m3) 

e 

n 

w 

(%) 

0 

(%) 

s 

(%) 

1/BD1 

127  2.03 

104 

1.66 

0.63 

0.39 

22.7 

37.6 

97 

2/BD2 

126  2.01 

103 

1.64 

0.65 

0.39 

22.3 

36.4 

92 

3/BD3 

126  2.02 

104 

1.67 

0.62 

0.38 

21.2 

35.2 

92 

Y  =  moist  unit  weight  (density) 
Yd  =  dry  unit  weight  (density) 
e  =  void  ratio 
s  =  saturation 


n  =  porosity 

w  =  gravimetric  water  content 
0  =  volumetric  water  content 


by  embedding  two  stacked  2-  x  8-in.  (5-  x  20-cm) 
planks  in  the  soil  (Fig.  10).  We  dug  a  trench, 
placed  the  planks  in  the  trench,  backfilled  around 
them,  tamped  them  down  and  removed  them 
before  the  first  freeze. 


Instrumentation  and  measurements 

Soil  freeze-thaw 

I  froze  the  soil  by  placing  three  refrigeration 
panels  on  the  bin  top  and  circulating  the  refriger¬ 
ant  through  them  (Fig.  11).  The  panels  rested 


a .  Digging  trench  for  planks. 


b.  Tamping  planks  into  trench. 


c.  Simulated  rill  before  FI. 


Figure  10.  Preparation  of  simulated  rill 
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Figure  11.  Soil  bin ,  during  a  freeze,  covered  with  pan - 
els  and  an  insulating  blanket. 


creases  dramatically  (Atkins  1989).  This  marked 
change  in  resistance  is  a  more  reliable  indication 
of  the  phase  of  soil  water  than  soil  temperature 
alone,  because  a  soil  can  be  0°C  when  the  water 
within  its  voids  is  still  unfrozen  (Atkins  1979). 

When  most  of  the  soil  water  between  a  ring 
pair  is  frozen,  the  resistance  recorded  by  that  pair 
tends  to  level  off  until  the  soil  begins  to  thaw, 
then  the  resistance  drops  just  as  rapidly  as  when 
it  froze  (Fig.  14).  When  thaw  is  complete,  the 
resistance  returns  to  a  before-frost  value.  Thus,  I 
determined  the  frost  penetration  rate  in  the  bin 
from  the  resistance  curves  that  showed  the  time 
when  freezing  began  and  when  the  frostline 
reached  the  depths  of  the  successive  ring  pairs, 
which  were  at  known  depths. 


Table  3.  Measurements  and  instrumentation. 


Condition 

Measurement 

Instruments 

Measurement 

accuracies 

Measurement 

frequency 

0°C  isotherm 

Temperature  (°C) 

RG  Thermocouples 

±0.1°C* 

Hourly 

Frost  depth 

Soil  resistance 
(millivolts) 

Resistivity  gages 

Absolute  units 
not  measured; 
change  is 
important 

Hourly 

Soil  surface 
geometry 

Distance  from  datum 
to  soil  surface  (mm) 

Millimeter  stick 

+1  mm  (est.) 

Before  and  after 
freeze  and 
during  thaw 

Soil  water 

Temperature-corrected 
dielectric  constant- 
converted  to  %  bv  volume 

Hydra  Probes+ 

±  0.015  to 
0.020** 

±  0.003++ 

Hourly 

*  Within  temperature  range  of  the  experiment,  -12  to  20°C. 
t  Vitel,  Inc.  (1994). 

**  Absolute  accuracy  with  no  soil-specific  calibration  performed, 
tt  Relative  accuracy  in  the  same  soil. 


about  5  in.  (13  cm)  above  the  soil  surface  and 
cooled  the  air  between  them  and  the  soil.  The 
freeze  (F)  portion  of  a  FT  cycle  is  the  period  when 
the  bin  air  temperature  (Ta)  was  0°C  or  below. 
Thaw  (T)  began  when  the  panels  were  removed 
and  the  bin  soil  was  exposed  to  the  ambient  air; 
thaw  ended  when  soil  resistances  showed  that  all 
frozen  soil  had  melted. 

Frost  depth  was  measured  by  two  Resistivity 
Gages  (RG)  (Table  3)  constructed  and  installed 
according  to  Atkins'  specifications  (1979,  1990). 
Each  gage  consists  of  a  PVC  rod  with  18  copper 
rings  (Fig.  12)  to  measure  soil  resistance  and  10 
thermocouples  to  measure  soil  temperature.  The 
RGs  were  mounted  to  the  bottom  of  the  bin  (Fig. 
12b);  the  final  soil  surface  was  even  with  the  RG 
tops  (Fig.  13).  The  RG  measures  soil  resistance 
and,  as  the  water  in  the  soil  between  two  rings  on 
the  rod  freezes,  the  resistance  of  that  soil  in- 


Four  thermocouples  measured  the  tempera¬ 
tures  of  the  bottom  of  the  panels  (Tp),  the  air  (Ta) 
between  the  panels  and  the  soil  surface,  the  soil 
surface  next  to  RG  C  (Tssc),  and  the  soil  surface  in 
the  middle  of  the  rill  (Tssr)  (Fig.  13).  These  ther¬ 
mocouples  and  those  in  the  RGs  were  from  the 
same  manufactured  wire  lot,  and  four  from  that 
lot  were  calibrated  in  an  ice  bath  to  determine 
their  absolute  accuracy.  Three  sets  of  measure¬ 
ments  were  taken  at  5-minute  intervals  while  the 
bath  was  stirred.  The  12  readings  averaged 
-0.035°C,  with  a  standard  deviation  of  0.032°C, 
and  maximum  and  minimum  readings  were 
0.027°C  and  -0.080°C  respectively.  Thus,  all  the 
thermocouple  temperatures  have  an  absolute  ac¬ 
curacy  of  ±  0.1°C  and  a  relative  accuracy  of 
0.01°C.* 


*  Personal  communication  with  K.  Knuth,  CRREL,  1995. 
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Water  content 

I  monitored  soil  water  redistribution  during 
each  FT  cycle  using  six  Hydra  probes  (Fig.  13). 
The  probes  measure  a  soil's  high  frequency  (50 
MHz)  complex  dielectric  constant,  which  is  made 
up  of  a  capacitive  and  a  conductive  electrical  re¬ 
sponse  (Vitel,  Inc.  1994).  The  capacitive  part  is 
most  indicative  of  soil  water  content,  while  the 
conductive  portion  reflects  predominantly  soil 
salinity.  As  a  soil  gets  wetter,  the  capacitive  re¬ 
sponse  increases  and,  with  appropriate  calibra¬ 
tion,  the  dielectric  constant  measurement  is 
directly  related  to  soil  water  content.  The  probe 
also  measures  soil  temperature  with  a  calibrated 
thermistor  in  the  probe  head.  The  temperature  is 
used  to  remove  most  of  the  temperature  effects 
on  soil  dielectric  constant  when  those  data  are 
converted  to  water  content. 

When  we  installed  the  probes,  we  placed  soil 
around  the  tines  (Fig.  15)  to  ensure  good  soil-tine 
contact,  as  recommended  by  the  manufacturer. 
We  dug  two  holes  in  the  surface  of  each  soil  layer 
and  placed  one  probe  horizontally  in  each  hole. 
The  Hydra  probes  read  from  35-37%  volumetric 
soil  water  before  the  first  freeze,  which  compares 
well  to  that  determined  from  the  bulk  density 
cores  (Table  2).  The  RGs,  thermocouples,  and 
Hydra  probes  were  connected  to  Campbell  data¬ 
loggers  (Fig.  16)  for  data  acquisition,  storage,  and 
relay  through  the  FERF  computer  data  collection 
system  (Knuth  1989)  to  my  computer. 

Surface  geometry 

I  used  a  millimeter  stick  to  measure  the  vertical 
distances  from  an  aluminum  bar  to  the  soil  sur¬ 
face  before  and  after  each  freeze  and  thaw.  I 
aligned  the  stick  with  a  vertical  mark  on  the  bar 
before  each  reading  to  ensure  that  I  was  measur¬ 
ing  the  true  vertical  distance  to  the  soil  surface 
(Fig.  17).  As  a  check  of  accuracy,  I  measured  a  dis¬ 
tance  when  the  stick  was  5°  out  of  vertical.  That 
measurement  was  1  mm  too  long,  but  the  5°  was 
obviously  not  vertical,  so  it  was  easy  to  hold  the 
stick  less  than  5°  out  of  vertical  during  each  mea¬ 
surement.  Thus,  a  conservative  estimate  of  error 
for  each  measurement  was  +1  mm. 

The  bar  was  bolted  in  place  above  the  soil  at 
locations  1,  3, 5,  and  7  (Fig.  17b).  The  vertical  dis¬ 
tances  were  measured  from  stations  5  cm  apart 
horizontally  along  the  bar,  except  over  the  rill, 
where  the  stations  had  to  be  closer  together  to 
adequately  define  the  rill's  cross-sectional  shape. 
Frost  heave  after  each  freeze  was  determined 
from  the  vertical  distance  at  four  out-of-rill  sta¬ 


tions  and  two  in-rill  stations  per  location.  The  dis¬ 
tances  were  also  used  to  construct  rill  profiles 
from  which  the  rill  hydraulic  radius  was  deter¬ 
mined.  I  originally  planned  to  measure  all  seven 
profiles  (Fig.  17a)  with  an  automatic  acoustic  pro¬ 
filer.  However,  the  profiler  did  not  work  properly 
across  the  rill's  sloping  sides,  so  I  had  to  measure 
the  profiles  by  hand  and  only  did  four. 

RESULTS  AND  DISCUSSION 
FT  Cycle  1 

Frost  formation  and  thaw 

The  panel  temperature  (Tp)  was  lowered  to  -20 
to  -30°C  at  1400  hours  on  Julian  day  146  (146.14) 
for  about  11  hours  and  was  held  between  -8  to 
-10°C  (Fig.  18).  The  first  freeze  (FI)  began  at 
147.15  when  Ta  reached  0°C  and  ended  233  hours 
later  when  the  panels  were  removed  at  157.08. 
The  soil  surface  at  RG  C  reached  0°C  35  hours 
(149.02)  after  FI  began  (Fig.  18).  The  0°C  isotherm 
penetrated  to  a  depth  of  about  6  in.  (15  cm)  at  RG 
C  during  FI  (Fig.  19),  an  average  of  0.04  in./hr 
(0.10  cm/hr)  (Fig.  20a).  The  isotherm  penetrated 
to  about  5  in.  (13  cm)  at  RG  D  (Fig.  19),  an  average 
of  0.03  in./hr  (0.08  cm/hr)  at  RG  D  (Fig.  20b), 
while  the  frostline  penetration  rate  averaged  0.03 
in. /hr  at  both  gages  (Fig.  20a,b).  The  first  thaw 
(Tl)  was  complete  at  158.23  at  RG  C  and  at  158.15 
at  RG  D  (Fig.  14),  39  and  31  hours,  respectively, 
after  the  panels  were  removed. 

The  liquid  water  decreased  dramatically  in  the 
surface  soil  as  most  of  the  water  froze,  but  only  a 
small  decrease  occurred  at  depths  of  7  in.  (18  cm) 
and  greater  (Fig.  21)  as  soil  water  was  drawn  to 
the  freezing  front.  During  FI,  the  liquid  water  at 
these  depths  never  froze.  The  similar  changes  in 
the  distribution  of  soil  water  diagonally  across 
the  bin  (Fig.  13)  suggest  that  they  represent  the 
soil  water  movement  out  of  the  rill  throughout 
the  bin.  Differences  in  the  heave  between  the  out- 
of-rill  and  in-rill  stations  suggest  that  soil  water 
movement  under  the  rill  may  have  been  different. 

Rill  hydraulic  geometry 

The  soil  outside  the  rill  heaved  at  all  16  mea¬ 
surement  stations  an  average  of  1.5  cm  during  FI 
(Table  4),  and  the  difference  between  maximum 
and  minimum  heave  was  1.4  cm.  Soil  in  the  rill 
bottom  heaved  at  six  of  eight  stations  by  an  aver¬ 
age  of  0.4  cm  and  the  difference  between  maxi¬ 
mum  and  minimum  heave  was  1.2  cm,  which  is 
less  than  the  soil  outside  the  rill.  This  is  some- 
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Ring 

Pairs 

A: -  1-2 

B: .  2-3 

C: -  3-4 

D:  - .  4-5 

E:  -  5-6 

F: -  6-7 

<3: -  7-8 

H: .  8-9 

I:  -  9-10 

J:  -----lo-n 

K:  - 11-12 

L: - 12-13 

M: - 13-14 

N: . 14-15 

O: - 15-16 

P: . 16-17 

Q: - 17-18 


c.  Probe  burial.  d.  Probe  in  position . 

Figure  15.  Hydra  probe. 


Figure  16.  Instruments  connected  to 
FERF  data  collection  system. 
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NORTH  END 


t 

Rill 


a.  Profile  locations. 

Figure  17.  Soil 


b.  Measurements  along  profiles, 
surface  geometry  measurements. 


Table  4.  Frost  heave  (cm).  Out-of-rill  meas¬ 
urements  were  made  at  stations  105, 90, 30,  and 
15  on  profiles  1  (PI),  3, 5,  and  7.  In-rill  measure¬ 
ments  were  made  at  stations  65  and  55  on  the 
same  profiles. 


Out  of  Rill 
(n  =  16) 

In  Rill 
(n  =  8) 

Max. 

Min.  Avg. 

Max. 

Min. 

FI 

2.2 

0.8  1.5 

1.2 

0 

0.4 

(P3) 

(P7) 

(P3 

(PI, 7) 

F2 

3.1 

1.8  2.5 

2.2 

0.8 

1.4 

(P5) 

(PI) 

(P5) 

(PI) 

what  surprising  in  that  Tssr  was  usually  0.1  to 
0.6°C  lower  than  Tssc  (Fig.  18)  and  I  would  have 
expected  that  this  lower  temperature  would  have 
established  a  steeper  temperature  gradient  below 
the  rill,  causing  more  soil  water  movement  to  the 
freezing  zone  under  the  rill,  resulting  in  greater 
in-rill  heave. 

The  rill  profiles  measured  2  hours  before  FI 
and  30  minutes  after  T1  began  (Fig.  22a,b)  show 
that  the  frozen  rill  sidewalls  were  generally  near¬ 
er  to  vertical  than  before  they  were  frozen,  which 
suggests  that  the  soil  along  the  rill  wall  near  the 
crests  heaved  more  than  that  near  the  rill  bottom. 
Heat  loss  from  the  soil  near  the  rill  crest  may  be 
enhanced  because  of  the  two  surfaces  through 
which  the  heat  can  flow  to  the  air.  Consequently, 
more  soil  water  could  have  been  drawn  to  the  rill 
crest  and  frozen. 


The  degree  of  FT-induced  rill  side  slope 
slumping  would  be  directly  related  to  the  amount 
of  soil  water  that's  drawn  into  the  freezing  zone 
as  frost  penetrates  into  the  soil.  The  amount  of 
frost  heave  is  a  measure  of  the  amount  of  ice  in 
the  frozen  soil,  which  is  determined  by  the 
amount  of  additional  water  drawn  to  the  freezing 
soil,  which  in  turn  is  determined  by  the  tempera¬ 
ture  gradient  in  the  soil  and  the  amount  of  avail¬ 
able  soil  water.  Upon  thaw,  the  additional  soil 
water  tends  to  saturate  surface  soils,  making 
them  susceptible  to  soil  flows  and  slumps.  Ice 
lenses  in  the  surface  soil  and  scattered  needle  ice, 
0. 1-1.0  cm  long,  on  the  soil  were  visible  at  the  end 
of  FI  and  F2,  which  suggests  that  sufficient  water 
had  been  drawn  into  the  freezing  soil  to  saturate 
the  surface  soils  when  thaw  began. 

Water  was  added  to  the  surface  soil  at  the  end 
of  each  freeze.  Before  the  panels  were  removed, 
warm  refrigerant  was  circulated  through  them. 
Otherwise,  cold  refrigerant  warming  in  the 
panels  would  expand  and  rupture  them.  When 
the  panels  heat  up,  the  needle  ice  that  forms  on 
their  undersides  during  a  freeze  (Fig.  23a)  par¬ 
tially  melts,  and  slush  ice  and  meltwater  fall  onto 
the  soil  (Fig.  23b,  24a).  This  additional  water 
(I  did  not  measure  its  volume)  would  simulate  a 
storm  of  very  wet  snow  falling  on  a  frozen  soil 
surface  and  the  degree  of  saturation  of  the  sur¬ 
face  soil  as  it  starts  to  thaw  would  likely  be  com¬ 
parable  to  that  which  may  be  found  in  the  field 
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Julian  Days 

a.  FTC1. 

Figure  18.  Temperatures  of  the  panel  (Tp),  air  (Ta),  soil  surface  at  mid-rill 
(T sSr),  and  soil  surface  next  to  RG  C  (Tssc). 


18 


146  147  148  149  150 


151 


152  153  154  155  156  157  158  159 
Julian  Days 


a.  FTC1. 

Figure  19.  Subsurface  soil  temperatures  (top — RGC;  bottom — RGD). 
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Julian  Days 


Figure  21.  Volumetric  soil  water  (%)  measured  by  the  Hydra  probes  at  0600  each  day 


Vertical  Distance  (cm)  Ver,ical  Distance  (cm) 


very  early  during  spring  thaw  in  a  temperate  cli¬ 
mate. 

Some  of  the  scattered  standing  water  on  the 
soil  surface  at  the  beginning  of  T1  slowly  drained 
into  the  thawing  soil,  contributing  to  its  instabil¬ 
ity,  and  some  flowed  over  the  rill  crest,  transport¬ 
ing  sediment  into  the  rill.  In  some  locations,  the 
saturated  soil  adjacent  to,  but  outside  of,  the  rill 
flowed  over  the  rill  crest  into  it,  forming  micro¬ 
drainage  channels  about  1  mm  deep  along  the  rill 
sidewalls. 


As  thaw  continued,  the  saturated  sidewall  soil 
flowed  and  slid  down  the  rill  in  variably  sized 
soil  masses  (Fig.  24c)  along  its  entire  length,  al¬ 
though  especially  in  mid-rill  (Fig.  24b).  Tension 
cracks  along  the  rill  crests  (Fig.  24b)  define  the 
blocks  of  soil  that  had  begun  to  slide.  The  more 
dramatic  cross-sectional  changes  resulting  from 
the  rill  soil  failures  occurred  at  profiles  3  and  5 
(Fig.  22c  and  d)  in  the  middle  of  the  bin,  rather 
than  at  profiles  1  and  7  near  the  bin  walls  where 
freezing  was  less  intense. 
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These  changes  in  the  hydraulic  geometry  of 
the  rill  are  the  main  aspect  of  the  experiment  be¬ 
cause  they  affect  the  hydraulic  radius,  R,  of  the 
rill,  which  would  partially  determine  the  velocity 
and,  thus,  the  erosivity  of  subsequent  flows  in  the 


rill.  The  rill  R  before  FI  at  the  four  profile  loca¬ 
tions  averaged  about  5.5  cm  (Fig.  25,  Table  5). 
After  Tl,  it  changed  insignificantly  at  profiles  1 
and  7  and  decreased  by  about  33%  at  profiles  3 
and  5.  Such  a  reduction  in  R  would  cause  the  vel- 


(T 1 )  (T2)  (Tl)  (T2)  (Tl)  (T2)  (Tl)  (T2) 


Figure  25.  Hydraulic  radius ,  R,  after  Tl  and  T2. 


Table  5.  Rill  hydraulic  geometry  measurements. 


Profile 

Bottom 

wtdth 

(cm) 

Top 

width 

(cm) 

Depth 

(cm) 

P 

(cm) 

A 

(cm2) 

R 

(cm) 

1 

19.0 

19.9 

13.2 

46.7 

257.5 

5.5 

T1 

1 

11.8 

27.9 

12.4 

42.3 

246.1 

5.8 

T2 

1 

0 

44.6 

10.4 

52.4 

231.4 

4.4 

3 

18.1 

20.2 

12.7 

43.5 

242.7 

5.6 

T1 

3 

0 

35.7 

8.6 

41.2 

154.3 

3.7 

T2 

3 

0 

64.5 

7.8 

66.0 

250.8 

3.8 

5 

16.7 

20.2 

12.7 

42.6 

233.6 

5.5 

T1 

5 

0 

31.1 

9.2 

38.6 

143.3 

3.7 

T2 

5 

0 

54.7 

8.1 

58.2 

220.6 

3.8 

7 

17.3 

19.9 

12.7 

42.9 

235.4 

5.5 

T1 

7 

11.5 

24.8 

11.5 

38.3 

209.0 

5.5 

T2 

7 

0 

41.8 

10.4 

48.7 

216.5 

4.4 

ocity  of  a  bankfull  flow  to  decrease  by  24%  using 
Manning's  equation  (eq  2)  when  keeping  rill 
slope  (S)  and  roughness  (n)  constant.  Thus,  in  this 
hypothetical  example,  the  erosivity  of  such  a  flow 
in  that  rill  could  decrease  significantly  after  one 
FT  cycle. 

FT  Cycle  2 

Frost  formation  and  thaw 

The  Tp  was  lowered  and  maintained  between 
-7  and  -13°C  (Fig.  18).  The  second  freeze  (F2)  be¬ 
gan  at  159.13  when  Ta  reached  0°C;  Ta  remained 
between  -4  and  -6°C  until  174.08  (355  hours  later) 
when  I  removed  the  panels  to  begin  T2.  The  soil 
surface  at  RG  C  reached  0°C  at  159.18  (Fig.  18b) 
and  the  0°C  isotherm  penetrated  to  a  depth  of 
about  10  in.  (25  cm)  at  both  RGs,  an  average  of 
0.03  in. /hr  (0.08  cm/hr)  (Fig.  20c  and  d).  The 
frostline  penetration  rate  averaged  0.03  in. /hr  at 
both  RGs,  as  well  (Fig.  20c  and  d).  T2  was  com¬ 
plete  at  177.09  at  RG  C  and  at  176.16  at  RG  D  (Fig. 
14),  73  and  56  hours,  respectively,  after  the  panels 
were  removed.  The  soil  at  RG  C  froze  deeper  and 
more  completely  than  that  at  RG  D  and,  thus, 
took  longer  to  thaw. 

Rill  hydraulic  geometry 

F2  lasted  122  hours  longer  than  FI,  and  the  Ta 
during  F2  was  generally  1  to  2°C  lower,  causing 
more  water  to  migrate  to  the  freezing  front, 
which  resulted  in  more  liquid  water  in  the  surface 
soil  being  frozen  longer  during  F2  (Fig.  21).  All 
stations  in  and  out  of  the  rill  heaved  and  the 
amount  of  heave  was  greater  during  F2  (Table  4) 
because  of  this  additional  frozen  water.  The  out- 
of-rill  soil  heaved  more  than  in-rill  soil,  just  as 
during  FI.  The  temperature  of  the  soil  surface  at 
RG  C  was  generally  0.1  to  0.3°C  lower  than  that  in 


mid-rill  during  F2  (a  reversal  of  that  measured 
during  FI),  and  possibly  the  temperature  gradi¬ 
ent  in  the  out-of-rill  soil  was  greater. 

During  F2  the  in-rill  and  out-of-rill  soil  surfac¬ 
es  remained  unchanged  at  profiles  3  and  5  (Fig. 
22d  and  e)  and  there  were  minor  changes  at  pro¬ 
files  1  and  7.  After  6  hours  of  thawing,  however, 
substantial  side  slope  collapse  had  occurred  and 
the  rill  was  significantly  wider  and  shallower 
(compare  Fig.  24b  and  26).  Some  6  days  after  T2 
began,  the  rill  had  filled  in  significantly  (Fig.  22f) 
and  the  changes  at  profiles  3  and  5  were  greater 
than  after  T1  (Fig.  22d). 


Figure  26.  Rill  shape  6-7  hours  after  T2  began. 

The  R  at  profiles  1  and  7  decreased  24  and  20%, 
respectively,  after  T2  (Fig.  25),  while  the  R  at  pro¬ 
files  3  and  5  was  virtually  unchanged.  These  de¬ 
creases  in  rill  R  would  reduce  a  hypothetical 
bankfull  flow  velocity  by  17  and  14%,  respective¬ 
ly.  F2  lasted  long  enough  to  cause  the  near-wall 
parts  of  the  rill  to  freeze  sufficiently  that  there 
was  substantial  soil  collapse  there  during  thaw. 
Note  the  less  dramatic  shape  changes  at  profiles  3 
and  5  during  F2.  In  regard  to  FT-induced  changes 
in  soil,  numerous  others  (Chamberlain  1973, 
1981;  Akram  and  Kemper  1979;  Formanek  et  al. 
1984;  Yong  et  al.  1985;  Van  Klaveren  1987;  Benoit 
and  Voorhees  1990;  Othman  and  Benson  1993) 
have  noted  that  the  most  dramatic  change  occurs 
after  one  FT  cycle,  with  lesser  change  following 
subsequent  cycles. 

LESSONS  LEARNED 

Rings  1  and  2  on  the  RGs  responded  appropri¬ 
ately  as  a  pair,  but  the  others  responded  two  pairs 
at  a  time  rather  than  one  pair  at  a  time  (Fig.  14). 
Two  ring  pairs  should  not  respond  together  be¬ 
cause  the  freezing  front  penetrated  at  0.03  in. /hr 
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Figure  21 .  Soil  resistance  changes  in  silt.  (After  Atkins  1989.) 


(0.08  cm/hr),  much  less  than  the  0.5  in. /hr  (1.3 
cm/hr)  required  for  the  frost  to  penetrate  be¬ 
tween  two  ring  pairs  in  1  hour  (the  time  between 
samples).  Knuth*  reported  that  this  erroneous  re¬ 
sponse  was  likely  caused  by  a  capacitance  prob¬ 
lem  in  the  cable  leading  from  the  RGs.  I  used  the 
data  from  only  one  of  the  paired  ring  pairs,  which 
reduced  the  resolution  for  locating  the  frost  line 
from  0.5  to  1.0  in.  (1.3  to  2.5  cm)  for  these  tests. 
Future  tests  will  be  run  with  the  RGs  wired  differ¬ 
ently  to  avoid  this  problem. 

Atkins  (1989, 1990)  points  out  that  the  general 
shape  of  the  temperature  and  resistance  curves 
during  freeze  and  thaw  should  be  similar  (Fig. 
27).  However,  Figures  14  and  19  show  a  precipi¬ 
tous  drop  in  soil  resistance  and  an  instantaneous 
rise  in  temperature  at  the  very  beginning  of  T1 
and  T2.  Obviously,  ice  in  the  soil  does  not  thaw  as 
quickly  as  these  figures  suggest.  The  problem 


*  Personal  communication,  1995. 


was  that  the  meltwater  that  falls  onto  the  soil 
when  the  panels  are  heated  seeped  down  be¬ 
tween  the  RGs  and  the  soil  to  a  depth  of  4  to  6  in. 
(10  to  15  cm).  The  rings  and  thermocouples  then 
measure  the  resistance  and  temperature  of  this 
water.  Atkins  (1989,  1990)  discusses  the  problem 
of  water  seepage  around  an  RG  and  suggests 
burying  the  RG  so  that  its  top  is  6  in.  (15  cm) 
below  the  soil  surface.  However,  I  need  to  mea¬ 
sure  frost  penetration  down  from  the  soil  surface, 
so  in  future  tests  I  plan  to  cover  the  top  of  the  RGs 
with  plastic  before  removing  the  panels  to  keep 
the  meltwater  away  from  them. 

Freeze  at  RG  D  was  less  intense  than  at  RG  C 
(Fig.  20).  RG  C  is  located  near  two  walls.  Maxi¬ 
mum  heave  always  occurred  at  profiles  3  and  5  in 
the  middle  of  the  bin,  the  minimum  always  at  1 
and  7  nearer  the  walls,  and  more  soil  water  froze 
in  the  middle  of  the  bin  than  nearer  a  comer.  Heat 
penetration  into  the  bin  may  have  been  greater 
nearer  the  bin  comers  because  of  the  proximity  to 
the  two  walls.  Profile  measurements  during  fol- 
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Figure  28.  Soil  temperature  measured  by  Hydra  probe  3  at  0600 
each  day  when  soil  was  frozen. 


low-on  experiments  will  be  made  exclusively  in 
the  middle  portion  of  the  bin  soil. 

The  absolute  values  for  Hydra  probe  soil  water 
contents  were  inaccurate  once  a  freeze  began.  The 
probe  uses  the  soil  temperature  it  records  to  cal¬ 
culate  a  temperature-corrected  real  dielectric  con- 
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Figure  29.  Dielectric  constants  vs.  temperature- 
corrected  constants. 


stant  of  the  soil  (Vitel,  Inc.  1994).  This  constant 
gives  the  water  content.  However,  the  probes 
were  wired  so  that  they'd  switch  off  and  on  to  take 
readings,  but  the  switch  failed  and  the  probes  re¬ 
mained  on.  The  probes  then  heated  up,  which 
warmed  the  soil  around  them,  and  the  tempera¬ 
ture  they  recorded  was  incorrect.  Probe  3  in  the 
surface  soil  between  RGs  C  and  D  (Fig.  13)  mea¬ 
sured  6-8°C  above  freezing  (Fig.  28)  when  the  sur¬ 
face  soil  here  was  frozen.  However,  the  relative 
values  for  soil  water  are  likely  accurate  because 
the  real  dielectric  constants  vs.  the  temperature- 
corrected  constants  recorded  during  these  tests 
are  closely  correlated  (Fig.  29).  The  switch  prob¬ 
lem  in  the  Hydra  probe  circuits  will  be  corrected* 
before  the  next  experiments. 

CONCLUSIONS 

These  initial  tests  suggest  that  the  pattern  of  rill 
change  attributable  to  FT-induced  infilling  that  I 
hypothesized  may  be  appropriate  and  that  such 
infilling  can  dramatically  alter  the  hydraulic  ge¬ 
ometry  of  a  rill,  which  would  substantially 
affect  the  velocity  of  flows  within  it. 

The  patterns  of  frost  penetration  and  thaw  and 
soil  water  redistribution  suggest  that  the  soil  in 
the  bin  froze  and  thawed  in  a  natural  manner, 
which  mimics  that  in  the  field.  The  heave  of  the 
soil  during  freeze  and  its  settlement  and  slumping 
during  thaw  appeared  to  reasonably  simulate  the 
response  of  an  undisturbed  field  soil  to  FT  cycles. 
Thus,  this  type  of  soil  bin  with  side  and  bottom  rein¬ 
forcement  to  withstand  vehicle  loading  can  be  used 
in  future  experiments  on  vehicular  ruts. 


*  Personal  communication  with  K.  Knuth,  CRREL,  1995. 
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